skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Li, Zhi-Yun"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract We analyze the dense gas kinematics in two class 0/I protostellar cores, Per 30 and NGC 1333 IRAS 7, in the Perseus Molecular Cloud to determine whether their velocity structures are indicative of rotation. We examine the hyperfine structure of the N2H+J= 1–0 transition by combining 3″ (900 au) Atacama Large Millimeter/submillimeter Array measurements with 9″ (2700 au) measurements from the Green Bank Telescope. We use theCASA Feathermethod to combine these data in order to maximize our sensitivity across spatial scales. We fit the N2H+spectra to constrain the centroid velocity of the gas at each pixel and use these values to calculate the linear velocity gradient and specific angular momentum within apertures centered on each protostar with radii ranging from 5″ to 60″. Our results indicate that the velocity structure probed by the N2H+emission is likely not a result of core rotation. These findings are consistent with other studies in the literature that indicate rotation is often not evident on scales ≲1000 au. We instead suggest that the velocity structure we see is a result of torques caused by irregular density distributions in these protostellar systems. 
    more » « less
    Free, publicly-accessible full text available May 22, 2026
  2. Abstract The vertical settling of dust grains in a circumstellar disk, characterized by their scale height, is a pivotal process in the formation of planets. This study offers in-depth analysis and modeling of the radial scale height profile of dust grains in the HL Tau system, leveraging high-resolution polarization observations. We resolve the inner disk’s polarization, revealing a significant nearside–farside asymmetry, with the nearside being markedly brighter than the farside in polarized intensity. This asymmetry is attributed to a geometrically thick inner dust disk, suggesting a large aspect ratio ofH/R≥ 0.15, whereHis the dust scale height andRis the radius. The first ring at 20 au exhibits an azimuthal contrast, with polarization enhanced along the minor axis, indicating a moderately thick dust ring withH/R ≈ 0.1. The absence of the nearside–farside asymmetry at larger scales implies a thin dust layer, withH/R < 0.05. Taken together, these findings depict a disk with a turbulent inner region and a settled outer disk, requiring a variable turbulence model withαincreasing from 10−5at 100 au to 10−2.5at 20 au. This research sheds light on dust settling and turbulence levels within protoplanetary disks, providing valuable insights into the mechanisms of planet formation. 
    more » « less
    Free, publicly-accessible full text available August 14, 2026
  3. ABSTRACT We investigate the dynamics of dust concentration in actively accreting, substructured, non-ideal magnetohydrodynamic wind-launching discs using two-dimensional and three-dimensional (3D) simulations incorporating pressureless dust fluids of various grain sizes and their aerodynamic feedback on gas dynamics. Our results reveal that mm/cm-sized grains are preferentially concentrated within the inner 5–10 au of the disc, where the dust-to-gas surface density ratio (local metallicity Z) significantly exceeds the canonical 0.01, reaching values up to 0.25. This enhancement arises from the interplay of dust settling and complex gas flows in the meridional plane, including mid-plane accretion streams at early times, mid-plane expansion driven by magnetically braked surface accretion at later times, and vigorous meridional circulation in spontaneously formed gas rings. The resulting size-dependent dust distribution has a strong spatial variation, with large grains preferentially accumulating in dense rings, particularly in the inner disc, while being depleted in low-density gas gaps. In 3D, these rings and gaps are unstable to Rossby wave instability, generating arc-shaped vortices that stand out more prominently than their gas counterparts in the inner disc because of preferential dust concentration at small radii. The substantial local enhancement of the dust relative to the gas could promote planetesimal formation via streaming instability, potentially aided by the ‘azimuthal drift’ streaming instability that operates efficiently in accreting discs and a lower Toomre Q expected in younger discs. Our findings suggest that actively accreting young discs may provide favourable conditions for early planetesimal formation, which warrants further investigation. 
    more » « less
  4. Abstract We present new Atacama Large Millimeter/submillimeter Array observations that, for the first time, detect hydrogen and helium radio recombination lines from a protoplanetary disk. We imaged the Orion Nebula Cluster at 3.1 mm with a spectral setup that covered then= 42 → 41 transitions of hydrogen (H41α) and helium (He41α). The unprecedented sensitivity of these observations enables us to search for radio recombination lines toward the positions of ∼200 protoplanetary disks. We detect H41αfrom 17 disks, all of which are HST-identified “proplyds.” The detected H41αemission is spatially coincident with the locations of proplyd ionization fronts, indicating that proplyd H41αemission is produced by gas that has been photoevaporated off the disk and ionized by UV radiation from massive stars. We measure the fluxes and widths of the detected H41αlines and find line fluxes of ∼30–800 mJy km s−1and line widths of ∼30–90 km s−1. The derived line widths indicate that the broadening of proplyd H41αemission is dominated by outflowing gas motions associated with external photoevaporation. The derived line fluxes, when compared with measurements of 3.1 mm free–free flux, imply that the ionization fronts of H41α-detected proplyds have electron temperatures of ∼6000–11,000 K and electron densities of ∼106–107cm−3. Finally, we detect He41αtoward one H41α-detected source and find evidence that this system is helium-rich. Our study demonstrates that radio recombination lines are readily detectable in ionized photoevaporating disks, providing a new way to measure disk properties in clustered star-forming regions. 
    more » « less
    Free, publicly-accessible full text available April 9, 2026
  5. ABSTRACT Recent high angular resolution ALMA observations have revealed rich information about protoplanetary discs, including ubiquitous substructures and three-dimensional gas kinematics at different emission layers. One interpretation of these observations is embedded planets. Previous 3D planet–disc interaction studies are either based on viscous simulations or non-ideal magnetohydrodynamics (MHD) simulations with simple prescribed magnetic diffusivities. This study investigates the dynamics of gap formation in 3D non-ideal MHD discs using non-ideal MHD coefficients from the look-up table that is self-consistently calculated based on the thermochemical code. We find a concentration of the poloidal magnetic flux in the planet-opened gap (in agreement with previous work) and enhanced field-matter coupling due to gas depletion, which together enable efficient magnetic braking of the gap material, driving a fast accretion layer significantly displaced from the disc mid-plane. The fast accretion helps deplete the gap further and is expected to negatively impact the planet growth. It also affects the corotation torque by shrinking the region of horseshoe orbits on the trailing side of the planet. Together with the magnetically driven disc wind, the fast accretion layer generates a large, persistent meridional vortex in the gap, which breaks the mirror symmetry of gas kinematics between the top and bottom disc surfaces. Finally, by studying the kinematics at the emission surfaces, we discuss the implications of planets in realistic non-ideal MHD discs on kinematics observations. 
    more » « less
    Free, publicly-accessible full text available December 12, 2025
  6. Abstract Due to dust grain alignment with magnetic fields, dust polarization observations of far-infrared emission from cold molecular clouds are often used to trace magnetic fields, allowing a probe of the effects of magnetic fields on the star formation process. We present inferred magnetic field maps of the Pillars of Creation region within the larger M16 emission nebula, derived from dust polarization data in the 89 and 154μm continuum using the Stratospheric Observatory For Infrared Astronomy/High-resolution Airborne Wideband Camera. We derive magnetic field strength estimates using the Davis–Chandrasekhar–Fermi method. We compare the polarization and magnetic field strengths to column densities and dust continuum intensities across the region to build a coherent picture of the relationship between star-forming activity and magnetic fields in the region. The projected magnetic field strengths derived are in the range of ∼50–130μG, which is typical for clouds of similarn(H2), i.e., molecular hydrogen volume density on the order of 104–105cm−3. We conclude that star formation occurs in the finger tips when the magnetic fields are too weak to prevent radial collapse due to gravity but strong enough to oppose OB stellar radiation pressure, while in the base of the fingers the magnetic fields hinder mass accretion and consequently star formation. We also support an initial weak-field model (<50μG) with subsequent strengthening through realignment and compression, resulting in a dynamically important magnetic field. 
    more » « less
    Free, publicly-accessible full text available July 16, 2026
  7. Abstract Circumstellar disk dust polarization in the (sub)millimeter is, for the most part, not from dust grain alignment with magnetic fields but rather indicative of a combination of dust self-scattering with a yet unknown alignment mechanism that is consistent with mechanical alignment. While the observational evidence for scattering has been well established, that for mechanical alignment is less so. Circum-multiple dust structures in protostellar systems provide a unique environment to probe different polarization alignment mechanisms. We present ALMA Band 4 and Band 7 polarization observations toward the multiple young system L1448 IRS3B. The polarization in the two bands are consistent with each other, presenting multiple polarization morphologies. On the size scale of the inner envelope surrounding the circum-multiple disk, the polarization is consistent with magnetic field dust grain alignment. On the very small scale of compact circumstellar regions, we see polarization that is consistent with scattering around sourceaandc, which are likely the most optically thick components. Finally, we see polarization that is consistent with mechanical alignment of dust grains along the spiral dust structures, which would suggest that the dust is tracing the relative gas flow along the spiral arms. If the gas-flow dust grain alignment mechanism is dominant in these cases, disk dust polarization may provide a direct probe of the small-scale kinematics of the gas flow relative to the dust grains. 
    more » « less
    Free, publicly-accessible full text available May 21, 2026
  8. Abstract The dust emission polarization spectrum—how the polarization percentage changes with wavelength—serves as a probe of dust grain properties in star-forming regions. In this paper, we present 89–214μm polarization spectrum measurements obtained from SOFIA/HAWC+ for three star-forming clouds: OMC1, M17, and W3. We find that all three clouds have an overall decreasing polarization percentage with increasing wavelength (i.e., a “falling polarization spectrum”). We use SOFIA and Herschel data to create column density and temperature maps for each cloud. We fit for the slope of the polarization spectrum at each sky position in each cloud, and using the Pearsonrcoefficient, we probe each cloud for possible correlations of slope with column density and slope with temperature. We also create plots of slope versus column density and slope versus temperature for each cloud. For the case of OMC1, our results are consistent with those presented by J. Michail et al., who carried out a similar analysis for that cloud. Our plots of polarization spectrum slope versus column density reveal that for each cloud there exists a critical column density below which a falling polarization spectrum is not observed. For these more diffuse sight lines, the polarization spectrum is instead flat or slightly rising. This finding is consistent with a hypothesis presented 25 yr ago in a paper led by R. Hildebrand based on Kuiper Airborne Observatory data. This hypothesis is that regions shielded from near-IR radiation are required to produce a sharply falling polarization spectrum. 
    more » « less
    Free, publicly-accessible full text available March 4, 2026
  9. Abstract We present JWST-NIRCam narrowband, 4.05μm Brαimages of the Sgr C Hiiregion, located in the central molecular zone (CMZ) of the Galaxy. Unlike any Hiiregion in the solar vicinity, the Sgr C plasma is dominated by filamentary structure in both Brαand the radio continuum. Some bright filaments, which form a fractured arc with a radius of about 1.85 pc centered on the Sgr C star-forming molecular clump, likely trace ionization fronts. The brightest filaments form a “π-shaped” structure in the center of the Hiiregion. Fainter filaments radiate away from the surface of the Sgr C molecular cloud. The filaments are emitting optically thin free–free emission, as revealed by spectral index measurements from 1.28 GHz (MeerKAT) to 97 GHz (Atacama Large Millimeter/submillimeter Array). But, the negative in-band 1 to 2 GHz spectral index in the MeerKAT data alone reveals the presence of a nonthermal component across the entire Sgr C Hiiregion. We argue that the plasma flow in Sgr C is controlled by magnetic fields, which confine the plasma to ropelike filaments or sheets. This results in the measured nonthermal component of low-frequency radio emission plasma, as well as a plasmaβ(thermal pressure divided by magnetic pressure) below 1, even in the densest regions. We speculate that all mature Hiiregions in the CMZ, and galactic nuclei in general, evolve in a magnetically dominated, low plasmaβregime. 
    more » « less
    Free, publicly-accessible full text available April 2, 2026
  10. Abstract Asymmetric and narrow infalling structures, often called streamers, have been observed in several Class 0/I protostars, which is not expected in the classical star formation picture. Their origin and impact on the disk formation remain observationally unclear. By combining data from the James Clerk Maxwell Telescope (JCMT) and Atacama Large Millimeter/submillimeter Array (ALMA), we investigate the physical properties of the streamers and parental dense core in the Class 0 protostar, IRAS 16544–1604. Three prominent streamers associated to the disk with lengths between 2800 and 5800 au are identified on the northern side of the protostar in the C18O emission. Their mass and mass infalling rates are estimated to be in the range of (1–4) × 10−3Mand (1–5) × 10−8Myr−1, respectively. Infall signatures are also observed in the more diffuse extended protostellar envelope observed with the ALMA from the comparison to the infalling and rotating envelope model. The parental dense core detected by the JCMT observation has a mass of ∼0.5M, a subsonic to transonic turbulence of M  =  0.8–1.1, and a mass-to-flux ratio of 2–6. Our results show that the streamers in IRAS 16544–1604 only possess 2% of the entire dense core mass and contribute less than 10% of the mass infalling rate of the protostellar envelope. Therefore, the streamers in IRAS 16544–1604 play a minor role in the mass accretion process onto the disk, in contrast to those streamers observed in other sources and those formed in numerical simulations of collapsing dense cores with similar turbulence and magnetic field strengths. 
    more » « less
    Free, publicly-accessible full text available May 21, 2026